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Abstract 

We set up and analyze a random matrix model to study energy 
localization and its time behavior in two chaotically coupled systems. 
This investigation is prompted by a recent experimental and theoreti- 
cal study of Weaver and Lobkis on coupled elastomechanical systems. 
Our random matrix model properly describes the main features of 
the findings by Weaver and Lobkis. Due to its general character, our 
model is also applicable to similar systems in other areas of physics - 
for example, to chaotically coupled quantum dots. 

1 Introduction 

The statistical features of coupled systems have attracted considerable inter- 
est in many branches of physics. Random matrix theory (RMT) has been 
successfully used in many of those investigations. RMT was founded by 
Wigner pQ. It is a schematic model [2j in which the Hamiltonian, or, more 
generally, the wave operator of the system, is replaced by a random ma- 
trix. The necessary prerequisite is that the system be sufficiently "complex," 
implying that the matrix elements of the Hamiltonian, or wave operator, 
calculated in an arbitrary basis, behave like random numbers. It has been 
shown that the spectral fluctuations in numerous different systems, if mea- 
sured on the scale of the local mean-level spacing, are very well modeled by 
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RMT; see the reviews in Refs. |3 IH E]- Due to the connection with chaos, 
one frequently refers to those systems as quantum chaotic which show cor- 
relations of RMT type. Similarly, systems are often referred to as regular if 
they lack spectral correlations. 

We consider two coupled systems. We assume that either the two systems 
are chaotic before they are coupled or that the coupling itself introduces 
chaoticity if the separate systems are regular. This scenario is equivalent 
to the breaking of symmetries, if only two values of the quantum number 
belonging to that symmetry are taken into account. The statistical features 
crucially depend on the strength of the coupling measured on the scale of 
the local mean-level spacing. Many studies have been devoted to this issue 
of chaotically coupled systems or, equivalently, to symmetry breaking. We 
mention isospin breaking in nuclear physics [HI U\, symmetry breaking in 
molecular physics S\, symmetry breaking in resonating quartz crystals [Oj 
and coupled microwave billiards 10.. While these studies addressed the 
spectral correlation, several investigations in nuclear physics IT2*1 1131 ITlj 
focused on the statistics of the wave functions and related observables in the 
presence of symmetry breaking or similar effects. In all these cases, RMT 
approaches in the spirit of the Rosenzweig-Porter model ^H] were successful. 

Sometimes observables in the time domain such as spectral form factors 
are more appropriate than the eigenenergy correlations functions [21 EH IU E] • 
This is so, for example, in the case of the presently much discussed fidelity; 
see Refs. [To*! [T7\ IT%1 H§] and references therein. Another example is the study 
of the energy spread in chaotic systems [201 El] . In the context of coupled 
systems, the time evolution of wave packets was investigated in Ref. j22j . 

In the present contribution, we study energy localization in two coupled 
systems in the time domain. This problem was addressed in a recent work 
by Weaver and Lobkis [23J who measured the time dependence of the wave 
intensity distribution in two coupled reverberation rooms. To this end, these 
authors recorded the time response to an elastic excitation of two coupled 
aluminum cubes. Moreover, they investigated the same problem theoretically 
and they numerically calculated the response in coupled two-dimensional 
membranes. In our study, we se tup and analyze an RMT model, based on 
the approaches in Refs. fH>\ [7J Its general character makes our model 
useful for similar problems in different physics contexts. In particular, we 
expect that our RMT approach also applies to coupled quantum dots. 

The article is organized as follows. In Sec. [2]we sketch the work of Weaver 
and Lobkis [231 • I n Sec. El we set up the RMT model and analyze it ana- 
lytically and numerically. We compare our results to those of Weaver and 
Lobkis in Sec. EJ Discussion and conclusions are given in Sec. El 
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Figure 1: Results from the experiment by Weaver and Lobkis. The energy 
ratio between the different rooms is plotted versus time (in ms). We note 
the different scales on the y axes. Reprinted from Ref. |2*3*] . 

2 Experiment and numerical calculations 

As we aim at a comparison with their findings, we present the work of Weaver 
and Lobkis [23] in some detail. Thereby, we also introduce the notation and 
conventions. The system studied experimentally consists of two aluminum 
cubes coupled by a solid connection, manufactured out of a solid aluminum 
block. The corners of the cubes were removed to desymmetrize the structure. 
This was done to ensure "chaotic" motion. Elastomechanical wave modes 
were excited in one room, and the response was measured in the other room. 
In this way, 16 different curves of energy intensity versus time were recorded, 
each in a small region around a different frequency. The results show that the 
energy does not always spread equally over the two rooms. If the coupling 
is weak, then the wave intensity is higher in the room where the initial 
excitation was performed than in the other room, regardless of how long one 
waits. Hence, the energy ratio never approaches unity. This deviation from 
the equipartition of the energy in the two rooms is referred to as energy 
localization. The resulting data are shown in Fig. [TJ and as expected there 
is localization in the bins of larger mean-level spacing, but not in the bins of 
small mean-level spacing. We will discuss these results further in Sec. EJ 

Two numerical studies were performed on membranes with rough bound- 
aries [22] ■ The dynamics of the system was in both cases governed by dis- 
cretized wave equations in each of the rooms, and the coupling was realized in 
different ways in the two numerical calculations, to which we refer as Nl and 
N2 in the sequel. In the first one, Nl, the connection had the form of a win- 
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dow between the two membranes similar to the situation in the experimental 
setup. In the second numerical calculation, N2, the rooms were separated, 
but springs were attached to a few different sites in rooms 1 and 2, thereby 
coupling those sites. In Nl and N2, a nonvanishing initial condition was 
given to one site in one of the rooms, and the response was calculated at dif- 
ferent sites in the other room. The resulting time series were cosine-bell time 
windowed to focus on a specific instant in time. Then, the time series were 
Fourier transformed and integrated over a small region in frequency to accu- 
mulate data around a certain frequency. As in the experiment, 16 different 
curves of intensity versus time were obtained around different frequencies. At 
different frequencies, the systems have different effective couplings. There- 
fore, one expects the time behavior of the different curves to differ in the 
degree of localization, as well as in the way in which this asymptotic satu- 
ration value is reached. Due to the differences in the coupling mechanism, 
there will also be differences between the results of Nl and N2. 

Moreover, Weaver and Lobkis performed an analytical model study. The 
elastic wave equation for the state of the system U (t) , say, is of second order 
in time. To focus on the response in a narrow interval around a certain 
frequency Q, the ansatz U(t) = u{t) exp(-iQt) is made with the assumption 
that u(t) varies slowly with time. This leads to a first-order differential 
equation in time for u(t) which has the form 



where C is the wave operator of the original second-order equation. The 
energies E\(t) and B^if) in rooms 1 and 2 at time t are defined as the to- 
tal probability density of finding the system state u(t) in one of the states 
ipik, i = 1,2, which are good eigenstates in room i: 



Strictly speaking, Ei(t) is no energy. Nevertheless, we find this terminology 
introduced in Ref. appropriate and use it as well, because Ei(t) measures 
the degree of motion in room i. If no energy dissipates into the surrounding 
environment, the total energy E = Ei(t) + E 2 (t) is conserved - i.e., indepen- 
dent of time. In the sequel, it is always assumed that the system is excited in 
room 1, and the energy is measured in room 2. Analytical solutions for two 
coupled states are presented in Ref. [23J by employing different statistical 
assumptions. 




c + n 2 
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3 Random matrix model 



We set up the model in Sec. 13 . 1 1 The connection to the two-level form factor 
is established in Sec. 13 .2^ and a 2 x 2 version of the model is evaluated in 
Sec. 13.31 We discuss numerical simulations of the RMT model in Sec. 13.41 
Finally, we comment on chaotically coupled regular systems in Sec. 13.51 



3.1 Setup of the model 



Spectral correlations in elastomechanics have been shown to be well described 
by RMT [21]. This is also true in the case of symmetry breaking [Hj, which 
is of direct relevance for the present study. Thus, RMT is also likely to 
be capable of modeling the time behavior of elastomechanical systems. As 
the first-order equation (JTJ) has proved to be a good approximation to the 
experimental situation, we also base our model on this Schrodinger type of 
equation. Thus, it is more natural to replace (C + Q 2 )/2Q by the random 
matrix H than to replace C itself. It turns out that this is indeed the best 
choice. 

The appropriate RMT model is an extension of the one employed in 
Refs. [7] H3*j . The random matrix H modeling the operator (C + Q 2 )/2Q 
reads 

„ r h, o i r o v l . . 

H = n u + a T/t n • I 3 ) 



The two matrices Hi, % = 1,2, model the uncoupled rooms 1 and 2. They 
are real symmetric and have random entries. We draw them from (two 
independent) Gaussian orthogonal ensembles (GOE's). As the rooms in the 
experiments and the numerical calculations Nl and N2 were of the same 
size, the level densities were also the same. This can be adjusted in the 
RMT model by giving the matrices Hi the same statistical weights and the 
same dimension N, such that 2N is the total dimension of H. The matrix 
V also has random entries. The strength of the coupling is measured by the 
dimensionless parameter a. It is sufficient to always assume a > 0. The 
statistical weight of V is chosen such that the total H is in the GOE of 
2N x 2N matrices for a = 1. 

We write the eigenvalue equation for the total Hamiltonian H in the form 



H^ n = uj n m n 



n 



2N 



(4) 



The eigenvalues u n and eigenvectors are functions of the coupling param- 
eter a. It is convenient to introduce the notation 



^2n 



(5) 
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where ty in , i — 1, 2, is the projection of \l/ n onto the subspace i We emphasize 
that ty in , % = 1,2, are functions of the coupling parameter a. The eigenvalue 
equations for the Hamiltonians Hi are written as 



l ^2n' l p2n 



11 

n 



N 
N 



(6) 



where the eigenvalues uji n and eigenvectors ipi n are not functions of a. This 
difference between \E^ n and ipi n is an immediate consequence of the fact that 
H depends on a, while Hi and H 2 do not. We will also use the notation 
^2n — {O^^nV for the corresponding 2iV-dimensional vector with zeros in 
the first N components. 

At time t = 0, the system is excited in room one such that the state of 
the total system can be written as 



S 



s 




(7) 



where s in room one is not specified in detail. We refer to 5* as to the source. 
The time evolution of the source is then simply 

u(t) = T(t)S , where T(t) = exp (iHt) (8) 

is the time evolution operator. Using the eigenvectors ip2n, the energy in 
room 2 is thus given by 



N 



N 



n=l 



E 2 (t,a) = J]u(i)i =J2( T (t)S)^2n$ n (T(t)S) 



n=l 




l N 



T{t)S . 



(9) 



The block matrix only contains the unit matrix ljy for room 2. We write 
E 2 (t,a) instead of E 2 (t) to emphasize the dependence on the coupling pa- 
rameter a in the RMT model. We will study averages E 2 (t,a) over the 
ensemble of matrices introduced in Eq. (jHJ). Occasionally, we will also con- 
sider averages over the direction of the source. This average is denoted by 
angular brackets (•••). 

The total energy E is trivially conserved in the framework of our model. 
We have 



E = E 1 (t,a)+E 2 {t,a) = S-S = s-a. (10) 
Hence, we can always construct Exit, a), once E 2 (t,a) has been calculated. 
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As all correlations have to be measured [21 El IH E] on the local scale of 
the mean-level spacing D, we introduce an unfolded time r = Dt. This is 
also the scale on which a acts; we introduce the unfolded coupling parameter 
A = a/D. The energy on the unfolded scale is then given by 

£ 2 (t,A) = lim E~ 2 (r/D,D\). (11) 

N^oo 

Small values of a have a large effect on the correlation functions if the mean- 
level spacing D is small as well. 

3.2 Relation to the two- level form factor 

We now derive an estimate for the time evolution of E2(t,a) which should 
apply to strong coupling strength- i.e., to a parameter a which is large on the 
unfolded scale of the mean level spacing. If one assumes that the source S 
comprises excitations into all states ipin, one may average over the direction 
of the source. We find, from Eq. ©, 

(E 2 )(t,a) = Btr^ T\t) 

where tr^ 11 ) is the trace over the (11) sector of the whole matrix - i.e., over 
the upper left block. The constant B results from the average. Expanding 
the time evolution operator in terms of the eigenvectors, 

27V 

T(t) =^* n exp(zu; n *)*t , (13) 

n=l 

we arrive after a short calculation at 

(E 2 )(t, a ) = Bj2 (*im • *m) (^2n ■ ^2m) exp [i(uj n - LU rn )t] . (14) 

n,m 

The vectors *&i n and ^ 2n are, according to Eq. (j3J), the projections of the full 
eigenvector \l/ n onto the subspaces corresponding to the two rooms. They 
depend on the coupling parameter a. These vectors coincide with the eigen- 
vectors for the matrices Hi and H 2 only for a = 0. Hence, they are only 
orthogonal in this special case. For all values a > 0, the scalar products 
• \l/i n and \l/ 2n • ^ 2m are neither zero nor given by 5 nm . The ensemble 
average over the eigenenergies u n and over the eigenfunctions ty n decouples 
only for the parameter values a = and a = 1- i.e., if the system either falls 
into two GOE's or is represented by one GOE. We now estimate the ensemble 




In 



T(t) , 



(12) 
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average of the energy (E 2 )(t, a) by making the approximation that the aver- 
ages over eigenenergies and eigenfunctions decouple. For strong coupling this 
should yield a reasonable result, and we will use this approximation in that 
limit only. We average over the eigenfunctions. For each term in the double 
sum in Eq. ((T3J) the product of scalar products gives the same contribution 
which we can take out of the sums. Thus, we find 



(E 2 ){t, a) = -B^exp [j(u n - u m )t] , (15) 

n,m 

where we absorbed the contributions from the average over the scalar prod- 
ucts into the constant B. The average over the eigenvalues is yet to be 
performed. Luckily, the average in Eq. ()15j) is recognized as the Fourier 
transform of the two-level correlation function El- This yields, on the 
unfolded scale, 

e 2 (r, A) = B [5(r/2n) + K 2 (t/2tt, A)] , (16) 

where the function K 2 (t, A) is referred to as the two-level form factor. The 
term 5(r) is due to the diagonal contributions in the double sum of Eq. (fT3j) . 
We notice that the conventions used here require a rescaling of time with 
a factor of 2tt. The expression ([TE|) is an estimate for the energy e 2 (r,\), 
exclusively in terms of the Fourier-transformed two-level spectral correlation 
function for the transition from two GOE's to one GOE. Unfortunately, the 
two-level form factor is not known analytically for all values of A. The result 
for one GOE, corresponding to A — > oo, reads |2j 

f 2r-rln(2r + 1) for < r < 1, 

K 2 (t,oo) = < _ . Ot + 1 r ^ -. (17) 

v ' ' I 2 - r In ^ _ I for r > 1. v ' 

We notice that the function e 2 {t, oo) has the limit properties 

£ 2 (0,oo) = 0, 
lime 2 (r,oo) = B. (18) 

T— >00 

The second property will be used to fix the scale in comparison with the 
results of Ref. The estimate (fTB^l should apply to large coupling - i.e., 
to parameter values A ^> 1. 



3.3 Two-by-two model 

Quite often, one obtains surprisingly good information about an RMT model 
by restricting it to the smallest possible matrix dimension such that the non- 
trivial specific characteristics of the model are still present. This is success- 
ful in the case of the nearest-neighbor spacing distribution for the Gaussian 
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orthogonal, unitary and symplectic ensembles (GGOE, GUE and GSE, re- 
spectively) [3 IH E] , but also for crossover transitions . Here, we proceed 
analogously. We obtain a 2 x 2 RMT model by setting N = 1 in Eq. (JHJ). It 
turns out convenient to absorb the parameter a into the matrix element V 
such that 

" H x V 

V H 2 

and to readjust the probability density function P{H) for the ensemble av- 
erage accordingly: 

V2 



H 



(19) 



P(H) 



7rv7ro- 



exp 



V 2 

-Hi - Hi - 2 V - 



(20) 



As before, the GOE is recovered for a = 1. We introduce eigenvalue and 
angle coordinates 





V ' 




V 


H 2 





cos if — sin (p 
sin <p cos <p 



00! 

Cg> 2 



cos (p sin if 
— sin if cos if 



[2V 



which implies that the time evolution can be written in the same form with 
eigenvalues exp(iuit) and exp(iu2t). For the integration measure one has 



dHidH 2 dV = -|c<ji — tjj<i\dtjj\du] 2 df . 



(22) 



The source in this 2x2 model is simply given by S = (1, 0). Thus, collecting 
everything, we find, with Eq. Q, 



E 2 {t,a) 



V2 



2tt 



exp 



duj\ J du 2 \ui — tu 2 \ J dp 

-oo o 
(lu 1 + lu 2 ) 2 [ljJ l -UJ 2 ) 2 



2 sin 2if 
cos 2ip H — 

cr 



■ 2 o . 2 ^1 — ^2 

sin Zip sin 



t . 



(23) 



We introduce x = {u)\ + uj 2 )/\/2 and y = {u\ — uj 2 )/\/2 as new integration 
variables, perform the x and if integrals, and arrive at 

E 2 (t,a) 



2(a + l) 



DO 

a J dyy exp [-y 2 (l + a 2 )] 



[l (y 2 (l - a 2 )) - h(y 2 (l - a 2 ))] cos 2ayt . (24) 
9 



Here, Io and Ii denote the modified Bessel functions of zeroth and first order, 
respectively. 

Formula ()24)1 is the final result in the framework of our 2x2 model. 
In Ref. |2Hj, solutions for two-state models based on Eq. ((TJ) with different 
statistical assumptions were given. Although the general behavior is similar 
to our 2x2 RMT result (j22), the analytical forms of these solutions in 
Ref. [23] are different. 

For obvious reasons, the unfolding of formula (J2l|) is meaningless. Nev- 
ertheless, experience with 2x2 model for the spacing distribution tells that 
meaningful statements can be achieved if the transition parameter - i.e., the 
coupling strength a in our case - is interpreted properly. Here, we can do 
that in the following manner. For large times t, the function E 2 (t, a) becomes 
constant, because it reaches its saturation limit. The latter will depend on 
a. Thus, comparing the saturation limit for the 2x2 model with that of the 
2N x 2N RMT simulation or with those of the experiment and the numer- 
ical calculations Nl and N2 allows one, in principle, to interpret a on the 
unfolded scale. The saturation limit is easily obtained from Eq. (J23J). Due 
to the Riemann-Lebesgue lemma we have 



because the cosine function in the integrand oscillates so rapidly for large t, 
that the integral gives zero in the limit t — > oo. 

We notice that there is no equipartition of the energies for t — > oo at 
a — 1. This may seem a bit unexpected, because H is a 2 x 2 GOE matrix 
for that parameter value. It simply reflects the need to properly interpret 
the parameters, as just discussed. We study a two-state system and compare 
the probability for a transition between the two states with the probability 
of staying in one state. There is no reason why these should be equal. In 
the N — > oo limit, however, we expect equipartition for the full GOE. As 
the coupling parameter has to be measured on the scale of the local mean 
level spacing, the limit iV — > oo corresponds to the limit a — > oo in the 2x2 
version, and in that limit the right-hand side of Eq. (J25J) tends to 1/2 and 
therefore to equipartition. 

3.4 Numerical simulations 

The numerical simulations of the RMT model are performed as in Refs. 
ITS] . The unfolding of the results, however, is not done in the standard way. 
To compare with the numerical calculations of Weaver and Lobkis, it is more 
convenient to first unfold the level densities and then "refold" the spectra 




a 



(25) 
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Figure 2: Energy ratio E2/E1 versus (dimensionless) time. Data from Nl 
of Ref. [23], bin 16, plotted as dashed line. The result from Eq. (fTHjl for 
A — > 00, B set by property 2 of Eq. ()18)) and a rescaling of time to make the 
two curves fit for early times, is plotted as a solid line. 

of our model with the level densities of Ref. [23! • Thereby we ensure that 
the mean-level densities of the RMT model acquire a form given by the 
Weyl formula for the billiardlike system of Ref. [23] • The appropriate Weyl 
formula for the modal density in a square membrane with Dirichlet boundary 
conditions reads [23! 

(") = |£ - 5 . (*> 

depending on the frequency Q. Here, A is the area of the membrane and L is 
the side length of one room. Moreover, c is the wave velocity and is chosen 
as unity here. Using twice the area of a room (A = 2 x 198 2 ) and a value of L 
which takes the roughness of the boundaries into account (L = 2x3x 198), we 
aquire an expression for the level density of the entire system. A comparison 
with the Weyl formula of Ref. [2Z! shows that the terms from the extra edges 
introduced to model the disorder cancel. Hence, we use the formula for a 
square membrane, which is precisely Eq. (j2fij) . This Weyl formula is employed 
to refold our RMT model for both numerical calculations Nl and N2. In the 
latter case, the system is not really a billiard, due to the springs used for the 
coupling. Nevertheless, the Weyl formula should be a good approximation 
to the real level density. 

For every simulation, we generate random matrices of dimension 2N x 
2N = 100 x 100, unfold, refold, calculate ^(t, a), and average over the re- 
sults of 800 such simulations. As the functions E 2 (t,a) are not measured 
on the unfolded scale, but on that of the Weyl formula, we do not introduce 
the notation e 2 (r, A) in the present context. Moreover, we notice that the 
numerical calculations Nl and N2 depend on the mean frequency Q in the 
bin under consideration. Accordingly, we arrive at a two-parameter family 
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Figure 3: Comparison between the 2x2 RMT model and the low frequency 
bins of N2 of Ref. j2Sj- The energy ratio E2/E1 is plotted versus (dimension- 
less) time. 

of time-dependent curves E 2 (t,a,Q). We now associate each bin in the nu- 
merical calculation performed by Weaver and Lobkis with the corresponding 
frequency Q. This leaves us with a one-parameter family of time-dependent 
curves E2n(t, a), /j, — 1, ... ,16, for each bin labeled by 

3.5 Chaotically coupled regular systems 

In the RMT model defined by Eq. (J3J) and in its subsequent analysis, we 
always drew the matrices Hi, i = 1,2, from (two independent) GOE's, im- 
plying that we model the two rooms individually as chaotic systems — before 
the coupling is considered. One can also assume that the two rooms are reg- 
ular systems before they are coupled. In that case, the matrices Hi, i = 1,2, 
would be drawn from Poisson ensembles |3]- We also did such numerical 
simulations. The results are qualitatively the same if the coupling - i.e., the 
matrix V - introduces enough chaos. The main effect is an adjustment of 
the scales. What matters for the qualitative behavior is the chaoticity of the 
total system. 

4 Comparison with experiment and numeri- 
cal calculations 

We now compare the results of our RMT model with those of Ref. The 
quantity studied in Ref. j2H| is mostly the fraction of energy in room 2, or 
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time x1Q 4 



(a) Bin 2. (b) Bin 5. 




5 10 



time x1Q 4 

(c) Bin 16. 

Figure 4: Energy ratio EijE\ versus (dimensionless) time. Data from Nl 
of Ref. plotted as dashed lines and from refolded numerical simulations 
using the RMT model as solid lines. We notice the different E2/E1 scale in 
plot (a). 
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time x1Q 5 

(c) Bin 16. 

Figure 5: Energy ratio E2/E1 versus (dimensionless) time. Data from N2 
of Ref. plotted as dashed lines and from refolded numerical simulations 
using the RMT model as solid lines. We notice the different EijE\ scale in 
plot (a). 
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Table 1: Coupling parameters a for the 2x2 model used in Fig. 3(a) The 
asymptotic saturation E 2 jE values for N2 in bins 2, 3, and 4 are taken 
from Ref. |2Sj- The corresponding coupling parameters are calculated from 
Eq. 



bin 


E 2 /E (asymptotic) 


a 


2 


0.35 


2.33 


3 


0.2 


0.67 


4 


0.14 


0.39 



E 2 (t)/E. In our figures, we plot E 2 (t) / Ei(t) instead. We first consider the 
estimate, Eq. (|17p. involving the GOE form factor and the result (|T6*j) of the 
2x2 RMT model and compare with the numerical calculations Nl and N2. 
As formula (fTtjj) was derived under the assumption of strong coupling on the 
scale of the local mean-level spacing, it should apply to the high-frequency 
bins of Nl. Anticipating the later extraction of the coupling parameter, we 
already now mention that indeed A > 1 in that bin. In Fig. |21 we compare 
Eq. (|T6*|) with data from bin 16 (the highest frequency bin) of Nl. A good 
description is obtained, although no equipartition is reached. This implies 
that the effective coupling A is large, but not very large. We recall that 
formally A — > oo corresponds to the strongest coupling a = 1 - i.e., to one 
single system. Expansion of the form factor for short times r reveals a linear 
short time dependence of the ensemble-averaged energy e 2 (t, A). This is in 
agreement with the analytical discussion of Ref. • 

We turn to the 2x2 RMT model and compare to bins 2, 3, and 4 of 
N2. The values of the coupling parameter a are determined from Eq. ({25} 
and given in Table ^ together with the E 2 /E saturation values of N2. The 
2x2 RMT model curves and data from N2 are shown in Figs. |3(a")1 and [3(b)| 
respectively. As expected, smaller values of a correspond to a higher degree 
of localization - i.e., to a stronger deviation from equipartition. The general 
behavior of the results from the 2x2 model stays the same for large values 
of a. We notice that the time scale is different, because the data from N2 
were not unfolded. The similarity shown in Fig. |3] between the 2x2 RMT 
model and the Weaver-Lobkis results j2S] is remarkable. This parallels the 
success of 2 x 2 RMT models for the spacing distributions. 

The RMT simulations yield a one-parameter family of curves for each bin, 
to be compared with the numerical calculations Nl and N2. We use visual 
inspection to determine, for each bin, which curve fits best to the numerical 
calculations Nl and N2. Typical results for some of the bins are shown in 
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Table 2: Coupling constants A on the unfolded scale resulting from the RMT 
simulation, determined for some bins of the numerical calculations Nl and 
N2 with mean frequencies fl: see Figs. (J2J) and (JSJ). 



bin 


Q 


A for Nl 


A for N2 


2 


0.1562 


0.4750 


3.2316 


5 


0.4686 


2.0222 


0.7754 


9 


0.8851 


3.0586 


0.7502 


13 


1.3016 


3.9093 


0.8441 


16 


1.6140 


4.2427 


0.8899 



Figs. |3] and 03 In the low-frequency bins of N2 there is a discrepancy; see 
Fig. |5(a)[ The RMT simulation does not overshoot its saturation value as 
clearly as the data of Ref. [2B1- This very large overshoot within a short- 
time interval is, however, borne out in the 2x2 RMT model; see Fig. El 
By the visual fit we determine the coupling parameter a. We measure it 
on the scale of the local mean-level spacing. In Table we list the in 
this manner obtained coupling parameter A for Nl and N2 for the bins under 
consideration. The A values extracted for Nl go up with higher bins, because 
the saturation value comes closer to equipartition. This effect is also visible 
for N2, but not so pronounced, because the saturation value does not change 
much for higher bins. 

We find good qualitative agreement between the RMT simulation and 
Nl and N2. Importantly, fluctuations between different random matrices 
are quite large, and the RMT curves presented here are averages of E2 over 
800 random matrices. Among those some show very much closer similarity 
to the curves of Weaver and Lobkis [23J, and we believe that most of the 
discrepancy in the low-frequency bins of Nl are due to fluctuations resulting 
from the specific choice of the boundary and perhaps other parameters in 
the numerical calculations of Ref. |23j . The reason for this interpretation of 
the discrepancy is the peculiar form of the fluctuations in bin 2 of Nl; see 
Fig. |4(a)l Since these lower bins have a lower level density, they should also 
be more sensitive to this kind of fluctuations. In the intermediate-frequency 
bins, we find that the numerical calculations Nl and N2 are well reproduced 
by the RMT model. In the figures for the high-frequency bin, Figs. |4(e)| 
and |5(ejl the RMT results are seen to deviate downwards for large times. It 
is conceivable that the short-time behavior strongly depends on the specific 
realization of the coupling while the long-time behavior is universal. 

The result from the experiment is shown in Fig. ^ The qualitative fea- 
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tures are well described by the RMT model. It can be seen that the system 
localizes in the low-frequency bins, but does not do so in the high-frequency 
bins. In other words, the system approaches equipartition of the energy. 
The localization behavior depends on the strength of the coupling measured 
on the scale of the local mean-level spacing, and the mean-level density is 
significantly higher in the high-frequency bins. It seems that the saturation 
value has not yet been reached on the time scale visible in Fig. ^ A closer 
inspection shows an initial behavior similar to the one in the results from 
our numerical RMT calculations and a slight upwards trend of the energy 
ratio towards the end of the time window studied. According to Ref. (23j 
this means that, first, the expected asymptotic value of E 2 /Ei is reached in 
the middle of the time window studied and that, second, there is another un- 
known effect acting on longer time scales adding to the energy spread which 
is likely to be due to the coupling of the setup to the environment. Thus, we 
do not attempt to extract a quantitative estimate of the coupling parameter 
a for the experiment. 

5 Conclusion 

We have set up and analyzed an RMT model to describe the time behavior 
of coupled reverberation rooms. This system shows localization effects under 
certain conditions. Within our RMT model, we gave an estimate of the strong 
coupling behavior which involved the two-level GOE form factor. Moreover, 
we studied the 2x2 version of our RMT model analytically for arbitrary 
coupling strength and performed numerical simulations for the 2N x 2N 
version. 

From the comparison with the work of Weaver and Lobkis, we conclude 
that the RMT model yields a good qualitative description. Moreover, we 
find an interpretation of the localization effect by relating it to the universal 
features of RMT models for crossover transitions. 

Formally similar RMT models have been studied in connection with sym- 
metry breaking. Then, the parameter a measures the degree of symmetry 
breaking. This is so for isospin breaking in nuclear physics |I30EJEI]> sym- 
metry breaking in molecular physics fSJ, and symmetry breaking in resonating 
quartz crystals The experiment that comes closest to the present situa- 
tion is the study of spectral correlations in coupled microwave billiards JU] • 
In this work and in the experiment of Weaver and Lobkis > & n excitation 
initially in system one would stay there for all times if the coupling is zero. 
This is formally analogous to a conserved quantum number. 

In the RMT model the parameter measuring the size of the connection 



17 



has a most natural counterpart: the root-mean-square matrix element due 
to the coupling measured on the scale of the local mean-level spacing. Thus, 
there are two ways of making the effective, dimensionless coupling parameter 
A small and to thereby introduce localization: either the original coupling 
parameter a which always has a dimension is small or the mean-level spacing 
D is made large. 

It might be surprising that no equipartition of the energy is seen for all 
coupling strengths even if one waits very long. This would be the expectation 
if one compares the system with two water basins coupled by a channel. 
Suppose that initially the two basins are empty and the channel is closed by 
a gate. Now one of the basins is filled with water and the gate is opened 
at t — 0. Obviously, the water levels in the two basins will be equal after 
sufficiently long time. The speed with which this equipartition is reached 
simply depends on the cross section of the channel. 

In the present case of the two coupled acoustic rooms the situation is 
different, because the wave character of the excitations has to be taken into 
account. Thus, the crucial parameter entering is the size of the coupling 
connection - i.e., its geometrical width - compared to a typical wavelength. 
The first waves after the excitation which come from system 1 into system 2 
enter a silent territory and cause the first excitations there. The next waves 
coming from system 1, however, encounter these first excitations in system 2 
and interfere with them constructively or destructively. This process contin- 
ues and, of course, after being excited in system 2, waves also travel back into 
system 1. For smaller times, this complicated dynamics certainly depends 
strongly on the realization of the coupling. This is clearly so in the numerical 
calculation of Ref . |2B] • Nevertheless, the loneg-time behavior, in particular 
the saturation limit, shows universal characteristics, consistent with general 
features of quantum chaotic systems; see Ref. [3] • In the energy domain, this 
is borne out in the fact that the correlations on smaller energy scales are 
described by universal RMT features, while system-specific properties show 
up on larger scales, leading to deviations from the RMT prediction. To avoid 
confusion, we emphasize that these universal RMT features include those in 
the presence of symmetry breaking. By system-specific properties we mean, 
most importantly, the scales set by the shortest periodic orbits. 

It is worthwhile to realize that the chaoticity of the individual subsystems 
before the coupling is not crucial, if the coupling itself introduces enough 
chaos. We tested numerically that the behavior of the corresponding model 
for two regular subsystems coupled chaotically shows the same qualitative 
behavior. The saturation value is reached slightly faster, which simply means 
that time is rescaled. Finally, we mention that our model is not restricted 
to elastomechanics. It would also apply to coupled quantum dots and other 
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coupled systems. 
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